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ABSTRACT 


A set of computer programs has been developed to estimate the 
takeoff and landing maneuver of a given aircraft. The program is 
applicable to conventional, vectored lift and powered- lift concept 
aircraft. Portions of the program may also be used to evaluate the 
static performance of these t^es of aircraft. The aircraft is 
treated as a point mass confined to motion in a vertical plane, 
and rotational dynamics have been neglected. The required input is 
described and a sample case is presented. 
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COMPUTER PROGRAMS FOR ESTIMATING AIRCRAFT 


TAKEOFF AND LANDING PERFORMANCE 
By Jeff V. Bowles and Thomas L. Galloway 
Ames Research Center, Moffett Field, Calif. 94035 


INTRODUCTION 


A set of computer programs has been developed to estimate the 
takesff and landing maneuver of a given aircraft. The program is 
applicable to conventional, vectored lift and powered-lift concept 
aircraft. Portions of the program may also be used to evaluate the 
static performance of these type aircraft. The aircraft is treated 
as a point mass confined to motion in a vertical plane, and the 
rotational dynamics have been neglected. The takeoff subroutines 
are a modification and simplification of an unpublished program by 
V. R. Corsiglia of Ames Research Center. 

The user is required to provide two subroutines which compute 
the total force coefficients along and normal to the flight path, and 
determine various required engine characteristics. 

This report describes the various subroutines and the required 
input, the equations used, and the computational techniques involved. 
Also included is ?. listing of the total program and an example case. 



NOTATION 


Symbol 

Fortran Name 



a 

n 

S(9) 

- 

acceleration normal to flight path (rad/sec 2 ) 

a. 

T(6), S(8) 

- 

acceleration along flight path (ft/sec 2 ) 

C x 

CX 

- 

force coefficient along flight path 

Cy 

CY 

- 

force coefficient normal to flight path 

ENP 

ENP 

- 

number of engines 

g 


• 

acceleration due to gravity (ft/sec 2 ) 

h 

HABS, S (7 ) 


absolute altitude (ft) 

i w 

EYEW, EYEWNG 

- 

incidence of wing (deg) 

LF 

XLF 

- 

load factor 

q 


- 

dynamic pressure (lb/ft 2 ) 

R/C 

sell) , ROC, RTCL 

- 

rate of climb (ft/sec,fpm) 

S 

SW, SWING 

- 

reference wing area (ft 2 ) 

T 

THRUST 

- 

thrust, net or gross, per engine (lb) 

V 

T(4), S(4) , VEL 


aircraft velocity (ft/sec) 

V 

X 

S(10) 

“ 

horizontal component of V (ft/ sec) 

V 

y 

sen) 

- 

vertical component of V (ft/sec) 

v . 

VI 

- 

engine failure speed, EAS (knots) 

V R 

VR 

- 

rotation speed, EAS (knots) 

w 

W, WG 

- 

aircraft weight (lb) 

w f 

,F, WFUEL 

- 

fuel flow (lb/hr) 

a 

ALPHA 

- 

angle of attack (deg) 

y 

S(5), GAMMA 

• 

flight path angle (rad, deg) 

6 f 

DELFD 

- 

flap deflection (deg) 

6 

s 

DELSPL 

- 

spoiler deflection (deg) 


v 



0 THETAF 

v ANGLE 

p RHO 

y MU 

P BRK XMUBRK 


- pitch attitude fuselage angle (deg) 

- vectored thrust deflection angle (deg) 

- air density (slugs/ft 3 ) 

- rolling coefficient of friction 

- braking coefficient of friction 


vi 



SUBROUTINES 


Subroutine TAKOFF 

The subroutine TAKOFF simulates the takeoff maneuver of a given 
aircraft. The program is applicable to conventional, vectored lift 
and powered lift aircraft. The aircraft is treated as a point mass 
confined in motion to a vertical plane and the rotational dynamics 
are neglected. This simplification necessitates an estimation of all 
rotational rates involved. These rates are either input by the user 
or are approximated by a finite difference form. 

The takeoff maneuver is divided into four basic segments: ground 

roll and rotation, liftoff and initial segment climb, acceleration to 
final climb speed at a constant rate of climb, and finally, the pullup 
maneuver to establish the final climb speed. Provisions in the program 
are made for gear retraction, flap retraction, changing of vectored 
thrust angle and changes in the power setting. 

The ground roll is made at a specified power setting and flap 
deflection. When the rotation speed is reached, the aircraft is 
"rotated" by increasing the angle of attack linearly with time until 
liftoff occurs or the tail scrape angle is reached. If the latter 
occurs, the ground roll is continued with the fuselage angle equal to 
the tail scrape angle. 

The flight path control is obtained by monitoring three dynamic 
variables — acceleration along flight path, load factor and fuselage 
angle (pitch attitude) - The aircraft is not permitted to decelerate 
and the load factor and fuselage angle are restricted to values less 
than or equal to a specified value. If any of these conditions are 
violated, the angle of attack is reduced until all constraints are 


satisfied. 



Once a specified altitude is attained, called the maneuvering 
altitude, the aircraft is pitched down by a reduction in angle of 
attack until a specified rate of climb is obtained. The aircraft 
then accelerates at this rate of climb until the desired final climb 
speed is reached. 

When the final climb speed is attained, the pullup maneuver is 
executed in order to bring the aircraft to a zero rate of acceleration 
along the flight path. This maneuver is accomplished by increasing 
the angle of attack and pulling a load factor of 1.20, which will 
result in an increase in the rate of climb to a final value at the 
desired climb speed. It may be necessary to throttle the engines in 
order to maintain the desired constant climb speed subject to the 
fuselage angle restriction. 

Program Inputs 

The inputs to subroutine TAKOFF are through the argument list, 
input by NAMELIST and common blocks /UNIV/ and /AERO/. 

The call to TAKOFF is as follows: CALL TAKOFF (INPC, IDCN, 


WGROSS, SWING, XENG, VI, VR, VEND) 

where 


INPC - program control 

- 

1 

- input data loaded 


a 

2 

- program executed 


s 

3 

- data input and 
program executed 

IDCN - print control 

a 

1 

- no print out 


• 

9 

- print out 


WGROSS - aircraft gross weight (lbs) 


SWING 


reference wing area (sq ft) 



XENC 


- number of engines 
VI - critical engine failure speed (knots) 

VR - rotation speed (knots) 

VEND - final climb speed (knots) 

All speeds are indicated air speeds. 

There are two namelist inputs to TAKOFF, /NAM1/ and /NAM2/. 
The namelist /NAMi/ input variables are as follows: 


CDGEAR 

- 

drag increment due to gear 

DFLPDT 

- 

flap retraction rate (deg/sec) 

DTABS 

“ 

temperature increment above S9°F (°F) 

DTGR 

- 

time required to retract gear (sec) 

DTPDWN 

- 

throttling down rate (percent/sec) 

DTPUP 


throttling up rate (percent/sec) 

DTVECT 

- 

vectored thrust angle reduction rate (deg/ sec) 

EYEWNG 

- 

wing incidence angle (deg) 

HAPT 

- 

airport altitude (ft) 

HGR 

- 

altitude at which gear retraction is 
started (ft) 

HMAN 

- 

maneuvering altitude (ft) 

HI (AX 

- 

takeoff termination altitude (ft) 

I OUT 

- 

for engine out takeoff, set IOUT * 1 

UM 

- 

rolling coefficient of friction 

NPAGE 

- 

no. of lines print od per page 

PMARG 

- 

pullup speed margin 

RTCL 

- 

rate of climb during accelerate segment (fpm) 

THTFLY 

- 

maximum allowable fuselage angle while airborne (deg) 

THTSCP 

- 

tail scrape angle (deg) 



XLFMAX 


maximum allowable load factor 


DTFAIL - time required for one engine to fail (sec) 

The user may input all, some, or none of the above input variables. 
The default values of these input variables are listed below: 


CDGEAR 

= 

0.0, DFLPDT = 

3.0 deg/sec, DTABS = 0. 

DTGR 

S 

5.0 sec, DTPDWN 

= 5.0 percent/ sec, 

DTPUP 

s 

6.0 percent/sec. 

DTVECT = 10. deg/ sec. 

EYEWNG 

= 

1.0 deg, HAPT = 

0.0 ft, HGR = 25.0 ft 

HMAN 

= 

1000. ft, I OUT 

= 0, UM = 0.02, PMARG 

RTCL 

s 

750. fpm, THTFLY 

= 15.0 deg. 

THTSCP 

s 

10. deg, XLFMAX 

= 1.15, DTFAIL = 3.0 


If the default value cf CDGEAR is used, the program will calculate, 
based on an empir. cal fo-*-:„,ia, a value for the gear drag as a function 
of gross weight and wing area. 

The second set of namelist variables, /NAM2/, constitute the 
flap, throttle and vectored thrust schedules. These are tables that , 

manage the flap setting, power setting and vectored thrust angle as a 
function of the aircraft speed and altitude. These variables are ! 

f 

arrays of dimension 5. 

| 

flap deflection (deg) 
flap retraction altitude (ft) 

flap retraction speed (knots) I 

power setting j 


XDELFD(I) - 
XHFLAP(I) - 
XVFLAP(I) - 
XPONER(I) - 



> 


l 


j 

i 

I 

* 

f 


XHPWR(I) - 
XVPWR(I) - 
XNU(I) 
XHVECT(I)- 
XWECT em- 


power setting change altitude (ft) 
power setting change speed (knots) 
vectored thrust angle (deg) 
vectored thrust angle change altitude (ft) 
vectored thrust angle change speed (knots) 




All altitudes are absolute altitudes and all speeds a^e indicated 
air speeds. 

These schedules are constructed as follows: If the speed or 

altitude of the aircraft is equal to, say, XVFLAP(I) or XHFLAP(I), 
respectively, then thj flaps are retracted to the value XDEI.FD(I). 

The power setting and vectored thrust angle management work in a 
similar manner. The power setting may either be increased or decreased. 

The flap setting and vectored thrust angle setting can only be reduced 
with speed and altitude. The values of XDELFD(l), XNU(l) and XPOWER(l) are 
all for the ground run. The user is permitted four changes in flap, 
power, vectored thrust angle settings during the airborne portion of 
the takeoff. 

The default values for /NAM2/ are as follows: 

- 100. percent throttle throughout takeoff 

- 0. degrees vectored thrust 

- 15.0 degrees flaps during ground roll, 
retracted to 5.0 degrees at 250 ft 
altitude, retraction to 2.0 degrees at 
200 knots, and finally, complete 
retraction at 210 knots 


I 

i 

i 

j 

I 

4 

l 

3 
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Again, the user may choose to use all, some, or none of the above 
schedule values. No changes to any of these settings are allowed 
during the pullup maneuver. 

Program Output 

The program output consists of a time history of several aircraft 
and flight path parameters. See sample listing presented in Appendix B. 
The variable DIST is the flight path track distance along the ground. 

The variable TAS is the true airspeed and EAS is the indicated airspeed. 
DALPH/DT and DTHTA/DT are the time rate of change of the angle of attack 
v.nd fuselage angle respectively. In addition, the user may also obtain 
the following values through the common block /EXCHNG/: 

SROLL - distance to liftoff (ft) 

S35 - track distance to 35 ft (ft) 

V35 - speed (EAS) sc 55 ft altitude (knots) 

T5J - ground distance covered by aircraft at 

time of engine failure 

The variable T5J may be used in conjunction with the subroutine 
ROLL to calculate accelerate — stop distances. 

The program will terminate normally when the end speed is reached 
(VEND) or when the maximum specified altitude (HMAX) is attained. 
Abnormal termination will occur under several conditions: 

- flight path constraints cannot be met by 
further reduction in angle of attack 

- aircraft cannot accelerate at input rate 
of climb (RTCL) 

- aircraft altitude goes negative 



- ground track distance greater than 
10. nautical miles 

- ground run exceeds 90 seconds 

- elapsed time greater than 300 seconds 

For further definitions a-"l explanations refer to the listing of 
TAKOFF and supporting subroutines contained in Appendix C, and the 
example case presented in Appendix B. 

Subroutine LANDING 

The subroutine LANDING simulates the landing maneuver of a given 
aircraft. As in the TAKOFF subroutine, the airex'aft is treated ; ' a 
point mass confined to motion in a vertical plane and the rotational 
dynamics are neglected. The landing is divided into three distinct 
phases: the approach, the flare to touchdown, and the rollout. 

The landing approach is made at the thrust and angle of attack 
required for zero acceleration along and normal to the flight path. 
These values of thrust and angle of attack for the steady state 
approach may either be input by the user or calculated in subroutine 
ZERO. 

The flare is executed by inci easing the angle of attack at a 
specified rate (DADT) . The program iterates on the flare initiation 
altitude until the rate of sink at touchdown equals a specified value 
(SINKTO) . 

Flight path control during the fl*;. . obtained by monitoring 
the load factor, fuselage angle and pitch rate. These three variables 
are all restricted to values less than or equal to values specified 
by the user. If any of these conditions are violated, the angle of 


attack is reduced until all constraints are satisfied. Provisions are 


made for spoiler retraction and increases in power setting during the 
flare if desired. 

After the touchdown, the rollout is then executed, with pro- 
visions in the program for any pilot reaction or delay time, brake 
application, deployment of spoilers and thrust reversing. The landing 
run terminates when the aircraft comes to rest. 

Program Inputs 

The inputs to subroutine LANDNG are through the argument list, 
common blocks /UNIV/ and /AERO/ and namelist /NAM3/. 

'Hie call to LANDNG is as follows: CALL LANDNG (INPC, xWRITE, 

WGROSS, SWING, XENG, PWRSET, MODE, IZERO, IREV, NER) where 


INPC 

- program control - 1 

input data loaded 


= 2 

program executed 


* 3 

input data loaded 
and program executed 

IWRITE 

- print control = 1 

no print out 


* 2 

print out 


» 3 

iteration on HFLARE 
printed out 

WGROSS 

- aircraft gross weight 

(lbs) 

SWING 

- reference wing area (sq ft) 

XENG 

- number of engines 


PWRSET 

- power setting 


MODE 

- roll out option control (see below) 



I ZERO 


- approach values of thrust and angle of 
attack = 0 calculated by program 

= 1 input by user 

IREV - reverse thrust control = 0 no reverse 

thrust 

- 1 reverse thrust 

used 

NER - error indicator = 1 successful case 


The user has two options for the ground roll calculations. An 
average deceleration rate and delay time may be input and the ground 
roll distance determined by a single formula (MODE = 2). The second 
option is the time step integration (MODE = 1) . 

The namelist /NAM3/ input variables are as follows: 


ABAR 

CDGEAR 

DADT 

DTABS 

DTDTMX 

DTPENG 

FAAFTR 

HAPP 

HAPT 

HSPOIL 

I POWER 

I SPOIL 

PWRIDL 

PWRMAX 


- average deceleration (g's) 

- drag increment due to gear 

- time rate of change of angle of attack 
during flare (deg/sec) 

- temperature increment above 59 °F (°F) 

- maximum allowable pitch rate during flare 
(deg/ sec) 

- throttling rate (percent/ sec) 

- field length factor 

- approach screen height (ft) 

- altitude of airport (ft) 

- minimum spoiler retraction altitude (ft) 

- throttle modulation control ) 


- spoiler retraction control 

- power setting at idle thrust 

- maximum power setting 


see below 
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PWRMRG 


- power setting margin 

SINKTD - sink rare at touchdown (ft/sec) 

TDELAY - delay time after touchdown before any 
braking action taken (sec) 

THEMAX - maximum fuselage angle (deg) 

TBRK - time before brake application (sec) 

TFLP - time before flap retraction (sec) 

TOFF - time before engine shut down (sec) 

TREV - time before reverse thrust applied (sec) 

TSPL - time before spoilers deployed (sec) 

XDIST - desiied landing distance (ft) 

XMUBRK - braking coeffient of friction 

XLFMAX - maximum allowable load factor 

QUESS - see subroutine ZERO 

HFLARE - initial quess of flare al. tude (ft) 

STEP - step size for iteration on HFLARE 

in addition, the following are input through common block /LAND/ 

COMMON /LAND/ GAMAPP, VKAPP, ALPHMX 

- flight path angle, measured positive 
below horizon (deg) 

- true airspeed at approach (knots) 

- maximum allowable angle of attack during 
flare, set equal to THEMAX in subroutine 
ZERO (deg) 

.V* 

V' * 

^ - 

JS;.\ : 




- 10 - 



The user may choose to input the approach speed (VKAPP) or input 
the desired landing distance (XDIST) and the program will calculate the 
required approach speed. This calculated approach speed is based on a 
constant load factor flare (XLFMAX) and the input values of ABAR, 

XLFMAX, GAMAPP, SINKTD and TDELAY. To use this option, the user should 
set VKAPP = 0. before the call to LANDNG. 

The approach values of thrust per engine and angle of attack may 
be input (IZERO = 1) through common block /AERO/ or calculated by 
subroutine ZERO (IZERO = 0) . 

The quess for the flare initiation altitude can be input *he 
user or calculated by the program (use default value of HFLARE) . The 
variable STEP controls the step size in the iteration to find the proper 
flare height. For powered — lift or any other aircraft that experience 
negative ground effects (suck-down) , it is suggested that HFLARE and 
STEP be input as 40.0 ft and 0.9 respectively. 

For ISPOIL = 0, the spoiler deflection angle (DELSPL) remains 
constant throughout the flare, and for IPOWER * 0, the throttle setting 
(PWRSET) is held fixed during the flare. For those aircraft experi- 
encing negative ground effects, the spoilers may be retracted to obtain 
direct lift control (ISPOIL ■ 1) and the power setting advanced during 
the flare (IPOWER ■ 1). If the load factor during the flare should 
decrease, usually occuring when at or near maximum allowable angle of 
attack, the power setting will be advanced. Spoilers may be retracted 
when the aircraft is at an altitude of HSPOIL or when PWRMAX - PWRSET < 
PWRMRG. 

If the user chooses to use reverse thrust during the ground roll, 
he must supply a subroutine REVRSE, which returns the total force 
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coefficients along and normal to the ground roll velocity vector. It 
is invisioned that this subroutine would be similar in structure to 
subroutine ARODYN, using the same common blocks to transfer required 
variable values. When using reverse thrust, IREV = 1 must be input. 

During the ground roll, the user may select to apply brakes, 
retract flaps, deploy spoilers, shut engines down or use reverse 
tnrust, or any combination thereof. Following touchdown, there is a 
delay of TDELAY seconds during which time no braking action is taken. 
After this delay time interval has elapsed, the various braking 
techniques listed above may be employed, with additional time delays 
for each (TBRK, TFLP, TOFF, TREV, TSPL) . Spoilers are extended at a 
rate of 90 degrees per second, 90 degrees being fully deployed, and 
flaps are retracted to zero degrees deflection at a rate of 10.0 degrees/ 
sec. If no reverse thrust is used, the power setting is reduced to 
idle (PWRIDL) at the specified rate of DTPENG percent/ sec. If reverse 
thrust is applied, the power setting is advanced to 100 percent at the 
specified rate of DTPENG percent/sec. If the user desires the flap 
settirg, spoiler deflection or throttle setting to remain fixed at the 
value at touchdown, he should input a large value for the associated 
time delay interval (say, 100 seconds) . 

The user may input all, some, or none of the input variables 
discussed above. The default values of these variables are listed below. 

ABAR ■ 0.35 g's, CDGEAR ■ 0.0, DADT * 5.0 deg/ sec, 

DTABS - 0.0 # F, DTDTMX - 7.0 deg/sec, DTPENG - 3.0 
percent/sec, FAAFTR ■ 0.60, HAPP ■ 35.0 ft, HAPT 
■ 0.0 ft., THEMAX « 10.0 deg, TBRK ■ 0.0 sec. 
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TFLP = 0.0 sec, TOFF = 0.0 sec, TREV = 0.0 sec, 



/ 

i 


i 

? 



*■ \ 



TSPL =0.0 sec, XDIST = 5000 ft, 

XMUBRK = 0.40, XLFMAX = 1.15, HFLARE = 0.0 ft, 

STEP =0.9 

If the default value of CDGEAR is used, the program will calculate 
a value for the gear drag based on an empirical formula as a function 
of gross weight and wing area. 

For further definitions and explanations refer to the listing 
of LANDNG and supporting subroutines contained in Appendix C, and the 
example case presented in Appendix B. 

Subroutine ZERO 

Given the aircraft configuration, its velocity and flight path 
angle, subroutine ZERO determines the required values of thrust per 
engine and angle of attack for zero acceleration along and normal to 
the flight path. 

The mathematical problem entails driving two functions of two 
independent variables simultaneously to zero. 

Let f (T,a) « dV/dt 
g(T,c.) « dy/dt 

where 

T - engine thrust 

a - angle of attack 

dV/dt - time rate of change of velocity along flight path 

dy/dt - time rate of change of flight path angle 


i 
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Subroutine ZERO generates the locus of f(T,ct) = 0 and g(T,a) ■ 0, 
and then searches for their intersection. 

The values of thrust and angle of attack at this point are the 
required value for the steady state approach. 



The above figure shows a typical plot of the f and g functions. 

The values a^ pp and T^ p are the desired values of angle of attack and 

thrust. The upper bound on the thrust is the static sea level thrust. 

The search for and is made on the interval ~ a max 5 5 ° max » 

i * 1,2. The error is then defined to be the difference between Oj and 

a . The subroutine varies the value of the thrust until the error is 
2 

driven to zero, or ■ a^. In varying the thrust, T^ + ^ ■ T^ times STEP, 
where STEP is an input parameter, usually less than 1.0. The initial 
quess at the value of thrust is defined to be the sea level static 
thrust times the parameter QUESS: 

T 1 - t sls #quess 
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It is suggested that the user make this initial quess of the 
thrust "high", so that the program does not search in areas where there 
is no solution of f(a,T) = 0. Recommended values of QUESS are: 

1) conventional aircraft — QUESS = 0.50 

2) powered — lift — QUESS = 0.80 

The default value for QUESS is 0.80. In the main calling program, the 
user must specify a value for the total sea level static thrust (ST) 
before the call to LANDNG is made. 

Subroutine ZERO may also be used to determine the static perform- 
ance of an aircraft, given its velocity and flight path angle. The 
flight path angle is measured positive below the horizon (i.e., if 
the aircraft is climbing at a flight path angle of 5.0 degrees, the 
call to ZERO would be made with GAMAPP = -5.0). 

The call to ZERO would be as follows: 


CALL ZERO (NER, EN, ALT, DTABS , KENG, PWRSET, QUESS) 

wher« 


NER 

EN 

ALT 

DTABS 

KENG 

PWRSET 

QUESS 


- error indicator * 1 successful case 

- aircraft Mach no. 

- aircraft altitude (ft) 

- temperature increment above 59*F ( # F) 
« 1 

- power setting (returned) 

- initial quess parameter 
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The following common blocks must also be used: 


/UNIV/ - passes ENP, ST, W, etc. 

/LAND/ - passes GAMAPP, ALPHMX 
/AREO/ - passes QS, VEL, HABS 

where 

QS = product of dynamic pressure and wing area 
VEL = aircraft speed (ft/sec) 

HABS = absolute altitude of aircraft (ft) 

Subroutines ENGINE and ARODYN are both called by ZERO. The call 
to ENGINE is made with KENG = 1 (see description of subroutine ENGINE) . 
The user may choose to work either with gross thrust or net thrust, 
as long as he is consistent in the use of the variable "THRUST" in the 
definition of total force coefficients and power setting (PWRSET) . 

Subroutines ARODYN and ENGINE 

The takeoff and landing subroutines described above require the 
user to provide two subroutines to compute total force coefficients 
and determine various required engine characteristics (e.g., thrust 
per engine and fuel flow) . The format and structure of these sub- 
routines is left to the discretion of the user. 

I. Subroutine ARODYN 

This subroutine computes the total force coefficients along the 
flight path and normal to the flight path as a function of angle of 
attack and thrust. A force coefficient in a particular direction is 
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defined to be the sum of all forces (excluding the weight component) 
in that particular direction divided by the dynamic pressure times the 
wing area. 


C 


x 


e 


- x 


l F 


qS 


The transfer of the various computer variables to and from sub- 
routine ARODYN is through labeled common blocks /UNIV/ and /AERO/. Of 
primary concern is the common block /AERO/: 

/AERO/ VEL, QS, HABS, THRUST, TVECT, ANGLE, DELFD, DELSPL, 
ALPHA, CX, CY, CL, CD, RHO, GRCD, I FAST 
The input variables from TAKOFF, LANDNG and ZERO are: 


VEL 

- 

aircraft velocity along flight path (ft/ sec) 

QS 

- 

dynamic pressure time wing area (lbs) 

HABS 

- 

absolute altitude of aircraft (ft) 

THRUST 

- 

thrust (gross or net) per engine (lbs) 

TVECT 

- 

total vectored thrust (lbs) 

ANGLE 

- 

angle of vectored thrust relative to aircraft 
center line, positive down (deg) 

DELFD 

- 

flap deflection (deg) 

DELSPL 

- 

spoiler deflection (deg) 

ALPHA 

- 

angle of attack (deg) 

RHO 

- 

air density (slugs/ft 3 ) 

GRCD 

- 

drag increment due to gear 


The return fron. ARODYN should be: 
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CX - total force coefficient along flight path 
CY - total force coefficient normal to flight path 

The output variable CL and CD are provided to the user as a means 
to distinguish between pure aerodynamic coefficients and total force 
coefficients. The output variables CL and CD are printed out in the 
time histories, but are not used in the actual calculations. If 
desired, in subroutine ARODYN, CL and CD may be directly equated to 
CY and CX, respectively. 

There is a certain amount of redundancy among some of the input 
variables. The user may utilize only those variables he desires and 
disregard the ethers. Due to the wide range of velocities encountered 
during the takeoff and landing maneuver, there will be a correspondingly 
large variation in the magnitude of the force coefficients which must 
be accommodated in subroutine ARODYN. 

II. Subroutine ENGINE 

This subroutine provides the various propulsion data to subroutines 
TAKOFF, ZERO and LANDNG. 

The inputs to subroutine ENGINE are through the argument list and 
labeled common blocks /AERO/ ruid /UNIV/. 

The call to ENGINE is as follows: 

CALL ENGINE (ALT, DTABS, EN, PWRSET, WFUEL, KENG) 

where 

ALT - aircraft altitude 

DTABS - temperature increment above S9*F (*F) 
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EN 

PWRSET 

WFUEL 

KENG 


aircraft Mach no. 
power setting (see below) 
fuel flow (lbs/hr) 

engine control parameter (see below) 


The variable PWRSET is defined to be: 


PWRSET = 


Net thrust 

Net thrust available 


and is the parameter used in controlling the thrust level. It is used 
for power setting management during the takeoff (e.g. throttle cutback 
for noise abatement), thrust modulation during the landing flare, and 
is calculated for the steady state approach. 

Used in conjunction with PWRSET is the variable KENG, the engine 
control parameter. For KENG = 0 or 2 (takeoff and landing respectively), 
the value of PWRSET is input and the thrust and fuel flow returned. In 
the static performance calculations of subroutine ZERO, the call to 
ENGINE is made with KENG ■ 1. Here, the thrust per engine is input and 
PWRSET and fuel flow are calculated. 

The user may choose to work with either the gross thrust per 
engine or the net thrust per -ngine, provided he uses the variable 
THRUST properly in the calculation of the total force coefficients. 

For example, when using gross thrust per engine, the ram drag must be 
included in the total summation of forces. If the gross thrust vector 
and ram drag vector are col linear, the user may choose instead to work 
simply with the net thrust. 

Refer to the example case presented for an illustration of 
subroutines ARODYN and ENGINE. 
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APPENDIX A 


EQUATIONS 


Equations used in TAKOFF and LANDNG 


1) Equation of motion along flight path 


dV/dt = (g/W) (-C x qS - W sin y) 


2) Equation of motion normal to flig. path 


dy/dt = (g/WV' CC y qS - W cos y) 


3) Equation of motion during ground roll 


where 


dV/dt - (g/W)[-Wy ♦ qS(C y vi - C x )] 


g « gravity constant 

W ■ aircraft weight 

q * dynamic pressure 

S n wing area 

flight path angle 
aircraft velocity 

total force coefficient along flight path 
total force coefficient normal to flight path 


y u 

V « 

c • 


4) Load factor 


qSCy 

XLF • a - -- 

W cos y 
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5) Constant rate of climb equation 


Rate of climb = ROC = V sin y 


For ROC to be constant with time. 


dROC 
dt = 


0 


dROC 

dt 


= ft: (V sin Y) 

= sin y + V cos y = 0 


Substitute for terms dV/dt and dy/dt from equations 1) and 2) 

(g/W)(-C x qS - W sin y) sin y 

+ V cos y ( g/WV) (C y qS - W cos > ) cos y s 0 

-C x qS sin y - W sin 2 y + CyqS cos y - W cos 2 y « 0 
or 

qS (C cos Y - C sin y) - W ■ 0 
Y * 

6) Rotational rate approximations by finite difference 


0 “ Y + o - i w 

where ( pitch attitude (fuselage angle) 

Y - flight path angle 
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m ipm urn * ** 



a - angle of attack 


i - incidence of wing 
w 


Differentiating with respect to time we obtain: 


d9 _ dot dy 

dt “ dt dt 


dv 

where is given by equation 2) . 

£ is approximated by the finite difference form: 


where 


da ( 
dt = la now 


a 

past 


)/At 


a 

now 

a 

past 

At 


= current value of angle of attack 
= previous value of angle of attack 
= integration time interval 


7) Simplified landing equations — for constant load factor (XLF max ) , 
constant speed (V^ pp ) flare 
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where: 


and 


R * 


V 2 

APP 


g(XLF 


max 


1 . 0 ) 


V APP 2 ( Y 


2 

APP 


- r 


TD 


flare 2g (XLF_ - 1.0) 


max 


XDIST * 


35.0 


R Y 


APP 


tan y 


APP 


\ y app / 


+ ' J>P t DELAY + 


'APP 


2g *a 


where 

t delay " delay time before braking 
a - average deceleration (g’s) 
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APPENDIX B 


SAMPLE CASE 


Shown below is an example of the input, calling format, sub- 
routines ARODYN and ENGINE, and the print out obtained for the takeoff 
and landing programs. 

The main calling program TEST1 is set up to do the takeoff and 
landing of a Boeing 727-200. The required common blocks are shown, 
but others may be added if required. Note that the spoiler deflection 
angle DELSPL has been set equal to zero before the call to TAKOFF. Also 
note that before calling the landing program, the flap deflection 
DELFD, the spoiler deflection DELSPL and the total static sea level 
thrust ST must be given values externally to LANDNG. 

The namelist input was as follows: 

$NAM1 NPAGE = 48, RTCL = 530., THTFLY = 20., HMAN = 2000. $END 

$NAM2 XPOWER{2) = 0.7S, XI!PKR(2) * 750., XPOWER(3) =0.95 
XHPWR(3) = 1750. $END 

$NAM3 DADT =1.0, IPOWER = 0, TFLP = 99., PWRIDL = 0.10, 

SINKTD = 3.0, TOFF « 3.0, TSPL = 2.00, XMUBRK = 0.45 $END 

Subroutine ARODYN calculates the lift and drag coefficients of 
the 727-200 as a function of angle of attack, flap and spoiler deflection. 
Th 3 increments of lift and drag due to flaps is determined in a table 
look-up format. The increase in drag and loss in lift increment due to 
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spoilers is assumed to vary linearly with spoiler deflection. Once the 
lift and drag coefficients are computed, the thrust components, normal- 
ized with respect to dynamic pressure times wing area (QS) , are adat. 
in to determine the total force coefficients CX and CY. 

Subroutine ENGINE computes the thrust and fuel flow of the JT8D 
engine, based on a simplified model. The thrust lapse is assumed to be 
linear with Mach number, and the fuel flow assumed linear with power 
setting. Note the use of the parameter PWRSET. The subroutine con- 
sists basically of two sections. For KENG = 0 or 2, the power setting 
PWRSET is input and the thrust and fuel flow returned. For KENG = 1, 
the thrust is input and PWRSET and fuel flow returned. No reverse 
thrust was used during the landing ground roll and r ibroutine REVRSE is 
simply a dummy subroutine. 

This particular run was made on the Lawrence Berkeley Laboratory 
CDC 7600 and required a CPU time of 2.38 seconds. 
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